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Accessing directly the properties of fundamental scalars in the confinement and Higgs 

phase 

Axel Maas^il 
Institute of Physics, Karl-Franzens- University Graz, Universitdtsplatz 5, A-8010 Graz, Austria 

(Dated: February 1, 2011) 

The properties of elementary particles are encoded in their respective propagators and interaction 
vertices. For a SU(2) gauge theory coupled to a doublet of fundamental complex scalars these 
propagators are determined in both the Higgs phase and the confinement phase and compared to 
the Yang-Mills case, using lattice gauge theory. Since the propagators are gauge- dependent, this is 
done in the Landau limit of the 't Hooft gauge, permitting to also determine the ghost propagator. 
It is found that neither the gauge boson nor the scalar differ qualitatively in the different cases. 
In particular, the gauge boson acquires a screening mass, and the scalar's screening mass is larger 
than the renormalized mass. Only the ghost propagator shows a significant change. Furthermore, 
indications are found that the consequences of the residual non-perturbative gauge freedom due to 
Gribov copies could be different in the confinement and the Higgs phase. 



I. INTRODUCTION 



Although the theory of weak isospin can be very well 
described with conventional perturbation theory after 
parametrizing the effective Higgs condensate [1|, it still 
poses quite a number of genuine non-perturbative ques- 
tions. The most obvious one is, what happens if the 
Higgs is very heavy. In this case, perturbative calcu- 
lations will break down eventually when reaching the 
energy domain of about 1 TeV [l|. But this is by far 
not the only question. A possibly even more interest- 
ing question is, whether the theory in itself can provide 
a dynamical mechanism for Higgs condensation, as in- 
dicated by lattice simulations [2|t4|- This would obliter- 
ate the necessity to drive electroweak symmetry breaking 
entirely by physics beyond the standard model, though 
not resolving the perceived hierarchy problem [8|. An- 
other serious obstacle is the possibility that the theory 
of weak isospin could indeed become trivial upon quan- 
tization, very much like the ungauged 0'* theory [9| or 
the ungauged Higgs- Yukawa theory [3, [l3|. However, 
the possibility of a non-perturbative stabilization (sim- 
ilar in concept, e. g., to asymptotic safety) has not been 
finally ruled out, and deserves further investigation. Of 
course, it could be that the top- Yukawa coupling will be 
important in case of the standard model 9] , but before 
requiring outside assistance, it should be better under- 
stood whether this is indeed necessary. Results concern- 
ing this question are not yet fully satisfactory. In par- 
ticular, besides the phase diagram [3, |^, Q usually onl 
gauge- invariant bound states have been investigated [2|, 
for this question, for which it may be difficult to separate 
genuine bound states from almost scattering states in a 
very weakly interacting theory. 



Aside from these questions, which concern the practi- 
cal viability of the weak isospin theory in its standard 
model form, there is also a more conceptual question. 
It has been shown that the confinement phase and the 
Higgs phase of the SU(2)-Higgs theory are not sepa- 
rated, if the theory is investigated with a (lattice-)cutoff, 
but are analytically connected [11|. Only in terms of 
gauge-dependent quantities can a phase transition be es- 
tablished, but the transition itself then turns out to be 
gauge-dependent 12|. This would imply that there is no 
qualitative difference between a confining and a Higgs 
phase in terms of gauge-invariant physics. Since a con- 
fining phase is not perturbatively accessible, so would 
neither be a Higgs phase. This must be sorted out to get 
a full understanding of the field theoretical setup of the 
weak isospin sector of the standard model. Furthermore, 
latest results indicate that the thermodynamic limit for 
questions concerning the phase diagram is rather hard to 
reach [6|]. 

On the other hand, this also permits to study an en- 
tirely different set of problems in the same theory. In 
the confinement phase, the color confinement process ap- 
pears to be entirely the same as for quarks, but without 
the complications introduced due to chiral symmetry and 
fermions^. For an understanding of color confinement, 
the weak isospin theory in the confinement phase there- 
fore offers an ideal laboratory, as has already been ex- 
ploited previously 14 18]. This is another incentive to 
study this theory in both phases. 

None of the points above, with the possible exception 
of a heavy Higgs, imply that for precision measurements 
in the standard model non-perturbative physics is neces- 
sary. Non-perturbative effects could easily be sufficiently 
suppressed over the whole accessible energy spectrum, 
even at LHC, or possibly up to the Planck mass, to make 
these questions practically irrelevant. And in fact, the re- 
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^ Note that the often cited Wilson criterion for confinement |13l| is 
actually blind to the spin structure of the involved particles. 



suits from LEP and Tevatron so far are in perfect agree- 
ment with a domination of perturbative contributions. 
However, as long as it is not clear what happens at LHC 
and beyond, the consequences of non-perturbative con- 
tributions should be known. 



To find answers to the questions above, in this work the 
possibility to access the properties of the (gauge-depen- 
dent) elementary degrees of freedom of the weak isospin 
theory will be investigated. These are the gauge bosons, 
which will be called the W for simplicity. Of course, since 
QED is neglected, the Z and W are degenerated here. 
The second will be denoted Higgs. Their propagators 
will be determined in both the Higgs and confinement 
phase, and will be compared to the quenched case, i. e. 
without dynamical Higgs particles, and thus a Yang- Mills 
theory with gauge group SU(2). For the non-perturbative 
determination a lattice implementation will be used. 



isospin modeP is given by [23 
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Since the propagators are gauge-dependent, it is neces- 
sary to fix a gauge. Given that non-perturbative gauge- 
fixing is made complicated by the Gribov-Singer ambi- 
guity [13, [201, a gauge is chosen in which these effects 
are comparatively well understood, which is the Landau- 
limit of the 't Hooft gauge, corresponding to the Landau 
gauge in Yang-Mills theory. Thus, as additional degrees 
of freedom, the Faddeev-Popov ghosts, also all degener- 
ate, are available. As a consequence of the Gribov-Singer 
ambiguity, on a finite lattice, as employed here, the Lan- 
dau gauge becomes a family of non-perturbative gauges 
in Yang-Mills theory [21] . This effect will also be studied 
here. 



This, rather exploratory, investigation is structured as 
follows: In section [H] the technical details of the simula- 
tions, including gauge-fixing, will be briefly discussed. In 
section lHll the propagators and their renormalization will 
be introduced and the central results presented. In sec- 
tion |IV] the consequences of the Gribov-Singer ambiguity 
will be analyzed. Everything is summarized in sectionFVl 
Some technical details are deferred to an appendix. 



In this expression a is the lattice spacing, W^ the gauge 
boson field, (j) the Higgs field, g the bare gauge coupling, 
A/k^ is the bare self-interaction coupling of the Higgs, 
Too the bare mass of the Higgs, the sums are over the 
lattice points x, and e^ are unit vectors on the lattice in 
the direction /i. 

For such a system the path integral, and thus the ex- 
pectation values of operators, can be evaluated using 
Monte Carlo simulations. The method used for the gauge 
part is described in detail in 24]. For the Higgs part, for 
each of the six sub-sweeps of the hybrid-overrelaxation 
cycle of the gauge field a local Metropolis update has 
been performed, where the acceptance probability was 
adjusted adaptively to be 50%. The results for observ- 
ables like the action, the spatial average of the expecta- 
tion value of the operator 



ri{x) = (f){x)^ (f){x) , 



(1) 



the expectation value of the plaquette, and the lowest 
masses of the isoscalar-scalar and isovector-vector exci- 
tations have been compared to literature values J2| to 
check the code. 

In the following, three different sets of parameters will 
be investigated, all for a lattice size of N* = 24^*^. One 
is the quenched case, i. e., the Higgs field is set to zero, 
while the other two cases correspond to systems deep 
inside the confinement phase and the Higgs phase'^ |15| . 



II. TECHNICALITIES 



The method used here is essentially standard lattice 
jauge theory. The lattice version of the weak SU(2) 



Note that strictly speaking the gauge group is SU(2)/Z2 due to 
the requirement of anomaly cancellation when including fermions 
[231 . Here, instead, the Z2 part is merely explicitly broken, and 
it is assumed that this does not make a difference on the level of 
correlation functions. 

Though these phases are not strictly separated on a finite lattice, 
it is possible to determine a phase diagram. Being sufficiently far 
away from the cross-over region, the systems are then classified 
by these names, though they are strictly speaking qualitatively 



A certain problem is posed by translating the lattice 
scale to a physical scale. In the Higgs phase it could be 
suggested to match with the Higgs condensate. However, 
this quantity is gauge- dependent, and in particular zero 
in the Landau gauge 12|. The possibility to instead use 
the spatially averaged expectation value of 77 is also not 
useful: It is non-zero also in the confinement phase 2], 
and it is therefore dubious to associate it with the Higgs 
condensate. The next possibility would be to match the 
W boson pole mass, which is measured in experiment. 
There are two reasons which make this approach rather 
hard. The first is that the pole positions would have to 
be determined by analytical continuation, which is rather 
unreliable with the number of lattice points available 
here. The second problem is that, at least in Yang-Mills 
theory, and likely also at least in the confinement phase, 
the W boson does not exhibit a pole |25l - [28| . Therefore, 
also this possibility is not useful for the present purpose. 
Concerning the Higgs, similar considerations apply, even 
if its mass would have been determined in experiment. 

Hence a more ad-hoc procedure will be used here. For 
both the confinement and the Higgs phase, the mass tti^ 
of the lowest lying state of the composite operator r]{x), 
^ , will be determined. This state would be an isoscalar- 
scalar Higgsonium bound state (or scattering state). Its 
mass will be arbitrarily set to 250 GeV, motivated that 
it may be twice as heavy as a single Higgs with its most 
possible mass of about 125 GeV 29]. Since, in contrast to 
quarks in QCD, it is found below that the Higgs screen- 
ing mass and its renormalized pole mass turn out to be 
rather close, it could be expected on the basis of a sim- 
ple constituent model that this should give an acceptable 
first guess. In the quenched case, this object is rather te- 
dious to calculate. Instead, here an upscaled version of 
the a(/3) relation of Yang-Mills theory will be used, to set 
a scale compatible with the one in the confinement phase. 
That said and done, one should see the physical scale to 
be rather of illustrative purpose, and it can always be 
scaled out again to replace it with another scale. 

It should be noted'* that in the lattice literature the 
operator ([T|) is usually associated with the Higgs itself 
rather than with a Higgs- Higgs-bound state 0, |j|; ^^^ 
thus its mass is denoted as the Higgs mass. The motiva- 
tion for this is that it is the simplest operator available, 
and, provided no anomalous hierarchy is encountered, 
it will also be the lightest state in the isoscalar sector, 
and thus the lightest physical excitation of the theory. 
In particular, it is the square of the radial mode of the 
elementary Higgs field, as being its lowest order gauge- 



identical. Furthermore, in the sense of the Wilson criterion, the 
confinement phase is not confining due to screening by pair cre- 
ation [1^, and in the Higgsphase the symmetry is actually not 
broken, but only hidden [1I. I12II . Therefore, these terms should be 
taken only as a short-hand notice, and their true meaning kept 
in mind all the time. 
* I am grateful to Christian Lang for a discussion of this issue. 



invariant contribution in an expansion in terms of poly- 
nomials in the fields. However, from the point of view of 
the gauge-dependent elementary fields, it is a composite 
bound-state, very much like the a-meson in QCD. It is 
thus not directly associated with the usual perturbative 
definition of the Higgs [1] in 't Hooft gauges where the 
Higgs is not an isoscalar. On the other hand, some dif- 
ferent relation may exist in the non-renormalizable uni- 
tary gauge. Anyway, it should be noted that neither the 
pole mass nor the screening mass of an elementary Higgs 
is a renormalization-group invariant, even if it should 
be gauge-invariant. In contrast, the mass of the low- 
est excitation in the rj channel is renormalization-group- 
invariant. Thus a direct identification of both concepts 
is, at least, not obvious. 

Part of the investigation here will therefore be to study 
how these two concepts are related^. One possible out- 
come could, e. g., be a similar relation as for quarks and 
mesons in QCD. 

That completes the basic idea used here to determine 
the scale. The details of the determination of the scale as 
implemented here are deferred to appendix |Al where also 
systematic uncertainties are discussed. The resulting set 
of parameters is then given in table HI The bare couplings 
have been chosen such that the systems are sufficiently 
far away from the phase transition such that no meta- 
stable phases could occur while at the same time the 
lattice spacing is not exceedingly small or large, based 
on the results in [2, 0, ll^l ■ 

This concludes the generation of the configurations. It 
remains to gauge-fix them. As stated in the introduction, 
this is potentially obstructed by the presence of Gribov 
copies. However, the presence of Gribov copies can also 
be turned into a virtue by using them to define different 
non-perturbative gau ges, depending on the selection of 
certain Gribov copies [2l|, [28[ [3^ . Still, all of these non- 
perturbative gauges satisfy the associated perturbative 
gauge condition, which therefore has to be chosen first. 

In the present case, this will be the Landau limit of 
the 't Hooft gauge. Thus, the gauge fields satisfy the 
condition 



df^W^ = 



(2) 



in the now Euclidean space time. To deal with the Gribov 
copies, first of all the selection of copies is restricted to the 
first Gribov horizon, defined to be the region with strictly 
positive semi-definite Faddeev-Popov operator [19| . 

The Gribov-Singer ambiguity now leads to the fact 
that potentially more than one gauge copy of a given con- 
figuration satisfies both conditions. In fact, in Yang-Mills 
theory appear to exist a large, possibly infinite, number 



^ Similar considerations apply to the lowest isovector vector chan- 
nel investigated in lattice calculations 0, Q , which is a collective 
excitation of elementary Higgs and W bosons, but the lowest 
gauge-invariant vector-particle state. 



TABLE I: The three sets included for comparison. For setting the scale a, see the text and appendix |X] The number of 
configurations are given for the determination of the propagators, where 1080 thermalization sweeps and 108 decorrelation 
sweeps have been performed in multiple independent runs. In the quenched case, the first number is for the W and ghost 
propagator, the second for the quenched Higgs propagator, mo is the tree-level mass of the Higgs. In the quenched case this is 
the mass appearing in the covariant Laplacian Q, see section Hill 
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of such copies in the continuum and infinite- volume limit 
[SlllS^]- At least in a finite volume, these can be used to 
design different gauges [2l|, |30| • In the next section, the 
so-called minimal Landau gauge [33| will be used, which 
selects one random representative among the remaining 
Gribov copies to represent a configuration. The conse- 
quences of alternative choices will be discussed in section 
HVl a detailed description can be found in [2l|, [23, l30ll. 

The method used to fix this gauge is described in [2J| . 
Since the gauge condition only involves the gauge fields, 
the same method for gauge- fixing as in Yang- Mills theory 
can be used. After obtaining the SU(2)-valued gauge- 
transformation g{x), the gauge- fixed Higgs field (j)' is ob- 
tained as 

(j)'{x) ^g{x)(p{x), 

as required to make the action gauge-invariant. 



III. PROPAGATORS 

Since Landau gauge is weak-isospin symmetric, and 
thus the vacuum expectation value of the Higgs field van- 
ishes ^^, the isospin symmetry is fully conserved and 
manifest even in the Higgs phase, as a consequence of 
Elitzur's theorem [3J|. Thus, also the propagators are 
isospin symmetric. There are then three independent 
ones, the W propagator, the Higgs propagator, and the 
Faddeev-Popov ghost [1[ . 

The W propagator in Landau gauge is given by 



7-jao ^ab 



P^ 



p^ 



The methods used to determine it and the ghost prop- 
agator below can be found in [2j|. In Landau gauge, 
one renormalization condition is required to fix the wave 
function renormalization constant^ -^w- It will be con- 
veniently chosen to satisfy 

,,2 

ZwZ{[?) 



^J■ 



M 



'w 



For simplicity, all wave-function renormalization constants are 
taken to multiply the propagators, not always in line with stan- 
dard conventions. However, since only the combination Z\yZ 
etc. play a role here, this is essentially irrelevant. 



with /i = 80 GeV and mw — 80 GeV. If the propaga- 
tor would be at tree-level, this would yield roughly the 
experimentally observed pole mass of the W boson. The 
resulting propagator Zy/Z(j?')l-p^ and the dressing func- 
tion ZwZ{p) are shown in figure [T] 

The results show that there is almost no difference be- 
tween the case without Higgs and the confinement phase. 
In particular, in both cases a non-zero screening mass ex- 
ists^, which is about l.lrriiY- For the Yang-Mills case, it 
is known that, despite its simple appearance, the ana- 
h^ic structure of the propagator is rather involved [25|- 
128( 1 . In particular, though a non-zero screening mass is 
present, there appears to be no pole mass, and the parti- 
cle has no representation as a Kallen-Lehmann state. A 
first glim pse indicates that the corresponding Schwinger 
function [23| in the confinement phase has essentially the 
same form as in Yang-Mills theory, suggesting that this is 
not changing. Nonetheless, this requires further investi- 
gation, in particular for larger lattice volumes and better 
discretizations. 

In the Higgs phase, the propagator is quantitatively 
different from the confinement phase. In particular, it is 
much closer to a tree-level behavior. Nonetheless, there 
are some differences, and in particular its screening mass 
is larger than the mass at the renormalization point, 
about l.lmw Therefore, there are sizable corrections 
to its behavior. Still, these are essentially quantitative 
effects^. Concerning the analytic structure, first indi- 
cations hint that at least quantitatively the Schwinger 
function in the Higgs phase may be substantially differ- 
ent. However, a qualitative difference to the confinement 
phase cannot be discerned unambiguously without larger 
and finer lattices. Therefore, this is not shown explicitly 
here. 

Instead of now directly investigating the Higgs propa- 
gator, it is interesting to investigate the ghost first. Its 



The actual value of the screening mass is very sensitive to lattice 
artifacts |35H37J |. and the value, which can be read from the plot, 
should be taken with care. 

Note that the apparent absence of a maximum of the dressing 
function in the Higgs phase is misleading. At sufficiently large 
momenta propagators in both phases will coincide, and both run 
logarithmically to zero essentially perturbatively [l||. Hence, in 
all cases there is a maximum, only its height and width change. 
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FIG. 1: The W propagator (right) and dressing function (left). Momenta are along the a;-axis for low momenta and along 
the a;y-diagonal for high momenta, such that violations of rotational symmetry for the displayed points are negligible. For 
larger momenta along the axes they become of the order of 10-20%. The quenched, confinement, and Higgs case are denoted 
by triangles, circles, and squares, respectively. Statistical errors in this section contain contributions from both statistical 
fluctuations and the scale uncertainty, the latter dominating. Note that the dressing function is dimensionless, and is therefore 
not affected by the scale uncertainty. 



propagator is a scalar function given by 






i^ab' 



P^ 



In the Landau gauge, its tree- level mass is zero \V\, and 
thus the renormalization condition 

ZaGi^i') = 1, 

will be used, with fj, — 80 GeV. The reason why this 
gauge degree of freedom is interesting is that in Landau 
gauge it is possible to combine it with the W propaga- 
tor to obtain an expression for the renormalization-group 
invariant running gauge coupling as [38|, |39[ 



a(/) = a{^l^)ZwZ{p^){ZGGip^)f 



(3) 



where the dependence on the renormalization scale of the 
propagators has been suppressed. This coupling is given 
in the so-called miniMOM scheme, but can be translated 
into the MS scheme, and is known up to four loops per- 
turbatively [40| . 

The result for the ghost propagator is shown in figure [5] 
It is immediately clear that there is a drastic difference 
between the confinement phase, again being essentially 
identical to Yang-Mills theory, and the Higgs phase. The 
significance of the infrared enhancement seen^ is likely 



associated with the confinement process. In marked con- 
trast is the ghost propagator in the Higgs phase, where it 
is almost a bare propagator, without any dressing. Indi- 
rect evidence for this drastic difference has already been 
obtained earlier from calculations in Coulomb gauge [44| , 
and has been discussed as a possibility in linear covariant 
gauges [451 - 147} . 

Hence, the most distinct difference between the con- 
finement and Higgs phase so far is the rather different 
ghost propagator. This also finds its manifestation in the 
effective coupling ([3]), as shown in figure [31 It is visible 
that the coupling has no Landau pole in the confinement 
phase, resembling the situation in the Yang-Mills case 
|35l437| once more. On the other hand, the coupling is in 
all cases infrared suppressed in this gauge, and starts to 
follow the same qualitative behavior at large momenta. 
Thus, from this point of view the only difference between 
the Higgs and the confinement phase is in the size of the 
coupling, rather than its low-momentum behavior. 

The final propagator is then the Higgs propagator. Its 
bare lattice version is obtained from the Fourier-trans- 
formed Higgs field 



i'b) = E' 



>'(a;). 



® At very small momenta the ghost dressing actually becomes also 
finite or logarithmically divergent |41| . However, the significance 
of this is currently still not fully resolved, see, e. g. |28l . |42| . |43| , 
for detailed discussions and in particular for further references. 



where p and x are lattice coordinates and momenta, and 
a symmetric lattice is assumed. The lattice Higgs prop- 
agator is then given by 
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FIG. 2: The ghost propagator (right) and dressing function (left). Momenta are along the a;-axis for low momenta and along 
the a;j/-diagonal for high momenta. The quenched, confinement, and Higgs case are denoted by triangles, circles, and squares, 
respectively. 
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FIG. 3: The effective coupling Q. Momenta are along various 
directions [2J|. The quenched, confinement, and Higgs case 
are denoted by triangles, circles, and squares, respectively. 
Note that the overall scale is arbitrary after division by q;(^). 



where the complex conjugation also inverts the direction 
of the momentum, as usual. This propagator can be as- 
signed physical units in the same way as for the W prop- 
agator to obtain the continuum propagator |24| . 

In the quenched case, the elementary Higgs field is 
not available. The quenched propagator is then obtained, 
in analogy to the quenched quark propagator, by the in- 
version of an operator. For a fundamental Higgs, this 
is just the fundamental covariant Laplacian with a mass 



term, 



D' 



da 



2 '^ 



"•O' 



with the generators of the gauge algebra t", which on 
the lattice is given by 



-Di = - 



M 



{u^,{x)&yi^.^+^^) + U+{x ~ ^i)5y^^,J,_^^) - 2i4yj 



+mQl5xy, 

where 1 is a unit matrix in weak isospin space. Since this 
operator is positive semi-definite, it can be inverted. For 
this purpose, the same method has been used as in case of 
the Faddeev-Popov operator in [2^. It should be noted 
that even a zero mass is not a problem for this method^'^ . 
With this, the bare scalar propagator is available for all 
systems. 

In contrast to the ghost and W propagator the 
renormalization of the scalar propagator is somewhat 
more complicated, since besides the multiplicative wave- 
function renormalization also an additive mass renor- 
malization is necessary. The renormalized propagator is 
given by [l| 



Dfip" 



sab 



Zh{p^ + vn?) + Hi:/(p2) + 5r 



^^ In contrast to the Faddeev-Popov operator, this operator has no 
trivial zero modes, and thus an inversion even at zero momentum 
is possible. However, since constant modes affect the result on a 
finite lattice, this is not done here. 
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FIG. 4: The Higgs propagator. Momenta are along the a;-axis 
for low momenta and along the a:j/-diagonal for high momenta. 
The quenched, confinement, and Higgs case are denoted by 
triangles, circles, and squares, respectively. 



phase, little is visible from the non-perturbative interac- 
tions. 

It should be noted that the propagators in the Higgs 
phase all exhibit the perturbatively expected behavior. 
In particular, a non-zero screening mass is found for the 
W and the negative mass squared of the Higgs became a 
positive effective mass. Thus, despite the vanishing vac- 
uum expectation value and thus the not explicitly hidden 
symmetry, the dynamics seem to be Higgs-like, and thus 
to be dynamically generated. However, one should be 
wary about this statement for two reasons: First, this is 
a result on a single lattice, and lattice artifacts may play 
a significant role. Secondly, this is not qualitatively dif- 
ferent from the confinement phase, again reemphasizing 
that the non-perturbative distinction of the Higgs and 
the confinement phase is still rather obscure. 

Finally, both the renormalized mass and the screening 
mass of the elementary Higgs agree with a very simple- 
minded constituent Higgs model for the Higgsonium state 
([1]). This therefore supports the view of the elementary 
Higgs building up the gauge-invariant states in the form 
of bound-states. 



where Hh is its self-energy, and Zh and Srri^ are the 
wave-function and mass renormalization constants, re- 
spectively. The two renormalization conditions imple- 
mented here are 
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where ^ = 125 GeV and the would-be pole mass mn is 
chosen to be 125 GeV. Using a derivative with respect 
to the momentum instead of the momentum squared 
is more appropriate for the lattice with its more-or-less 
evenly spaced momenta at the most reliable intermediate 
scales than the more conventional momentum squared. 
Of course, this is also sufficient to determine the two 
unknown renormalization constants, though the result 
is therefore a bit different from the usually employed 
schemes. 

The results for the renormalized Higgs propagator are 
shown in figured It is immediately visible that the prop- 
agator in all cases deviates only slightly from an almost 
bare propagator with the renormalized mass. The only 
significant deviations from the bare propagator are found 
at small momenta, and the screening mass Dh{0)^^'^ is 
in all cases larger than the renormalized mass. In the 
Higgs phase, its value is only slightly larger, about 130 
GeV, while in the confinement phase it is about 145 GeV. 
Thus, there remain some non-trivial infrared modifica- 
tions of the Higgs propagation. Other than that, the 
Higgs propagator is essentially unaffected compared to 
its tree-level behavior. Thus, even in the confinement 



IV. NON-PERTURBATIVE 
GAUGE-DEPENDENCE 

As noted, the Gribov-Singer ambiguity implies that in 
principle there exists more than one solution to the gauge 
condition ([2|) . The gauge used so far reduces the number 
of Gribov copies by fixing to the first Gribov region in 
Landau gauge, and then choosing randomly a represen- 
tative from the residual gauge orbit left. This implies 
that in this minimal Landau gauge the obtained gauge- 
dependent correlation functions are effectively averaged 
with a certain weight over the residual gauge orbits |2l| . 
To the best of our current knowledge, this weight func- 
tion is fiat, and thus the distribution on the gauge orbits 
is faithfully reproduced. 

This is only one admissible way to treat the residual 
gauge freedom. Many alternatives have been investi- 
gated, see e. g. [H [30, 48 50] . In Yang-Mills theory, at 
least at finite volumes and discretizations the correlation 
functions depend on this non-perturbative choice, in gen- 
eral up to the typical scale Aym ■ Though of course this 
dependency is irrelevant for the determination of gauge- 
invariant physical observables like cross-sections, it can 
be useful to choose an alternate gauge for technical rea- 
sons. E. g., there exists hints that some non-perturbative 
gauges in Yang-Mills theory may be more amendable for 
the explicit construction of the Hilbert space, while oth- 
ers, due to a lack of infrared singularities, are more useful 
in various practical calculations. It is therefore worth- 
while to study the dependence of correlation functions 
on this gauge choice also in the case with matter fields. 

For that purpose, the methods described in [2l|, |4g| 
will be adopted here. In particular, Gribov copies are 
then searched for by a multi-start algorithm in the gauge- 
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fixing procedure. In the present case, five random starts 
are performed for eacfi configuration, since only an in- 
dicative result is desired here. It is inherent to this 
method that only a lower limit to the gauge dependence 
of the propagators can be achieved. 

The first finding is that no Gribov copies are found 
in the Higgs phase. Besides the interesting option that 
indeed in the Higgs phase only configurations contribute 
to the path integral significantly which have no Gribov 
copies inside the first Gribov region, the second pos- 
sibility is that the strongly volume and discretization- 
dependent [21| number of Gribov copies for the se- 
lected parameters is just so small that none have been 
found. Given the experience with results on small two- 
dimensional volumes in Yang-Mills theory [2l|, [30| , where 
also very few copies are present, and the appreciable vari- 
ation of the propagators in the investigated momentum 
window here, either possibility would indicate that the 
residual gauge orbit in the Higgs phase has a significantly 
different structure than in the confinement phase. 

An alternative would be that the employed multi-start 
algorithm is just not successful in finding copies. A gain , 
given its successes for the Yang-Mills theory [2l|, |43, Ho] , 
this would imply a significantly different structure of the 
residual gauge orbit in the Higgs phase. 

Irrespective of the precise reason the residual gauge de- 
pendence of the correlation functions in the Higgs phase 
is thus not existent, and in the following only the con- 
finement phase will be investigated, and compared to the 
Yang-Mills case. In that case, the average number of 
Gribov copies found almost saturates the number of Gri- 
bov copies checked, i. e., there have been found 4.38(3) 
copies in the confinement phase per configuration in five 
attempts. In the Yang-Mills case, this number is slightly 
larger, 4.55(6). 

To provide an estimate of the variability of the correla- 
tion functions, these will be determined for the absolute 
Landau gauge and the max-i? gauge, see [21| for their 
definition. These gauges have been found so far to limit 
the variability of the correlation functions in Yang-Mills 
theory at a fixed number of Gribov copies. The first of 
these gauges attempts to minimize the W propagator by 
the choice of the Gribov copy, while the second attempts 
to maximize the ghost propagator in the infrared. 

The result for the three correlation functions and the 
running coupling are shown in figure [SJ The dependence 
of the quenched Higgs propagator is left out intention- 
ally, since it does not represent a dynamical variable, 
and therefore its gauge dependence is irrelevant for the 
present purpose. 

For the ghost propagator the gauge variability exceeds 
the statistical error. Thus, it is significantly gauge-de- 
pendent, even given the small numbers of Gribov copies 
used. In fact, its dependence is expected to increase sig- 
nificantly even further when enlarging the search space 
|50l [Slj . as is the case in Yang-Mills theory, if the ra- 
tio of search space to average number of Gribov copies 
found is so close to one as in the present case. On the 



other hand, the W propagator shows no change within 
the statistical error, though when including more Gri- 
bov copies this is known to change for pure Yang-Mills 
theory at fixed lattice parameters [48'|. Nonetheless, the 
running coupling, derived from the gauge propagators, 
is consequently significantly dependent on the gauge. In 
contrast, the Higgs propagator is not gauge-dependent 
within the statistical errors. This is not surprising, as 
the gauge conditions include only the ghost and W fields. 
However, if a 't Hooft gauge outside the Landau limit is 
used, this may change. Still, the gauge-dependence in 
the case with Higgs is rather similar to that in the Yang- 
Mills case for all correlation functions, again emphasizing 
the similarity between both cases. 

In total, the non-perturbative gauge-dependence in the 
confinement phase is sizable, and should be kept in mind 
if, e. g., results should be transferred between different 
calculations or even methods. On the other hand, within 
the limited scope of this investigation, no effect in the 
Higgs phase is found. 



V. SUMMARY 

Summarizing, the present study shows the feasibility 
of determining the gauge-dependent correlation functions 
in Yang-Mills-Higgs systems beyond perturbation theory, 
in both the confinement and the Higgs phase, even for a 
very heavy Higgs. It therefore provides a route to unify 
the non-perturbative gauge-invariant description usually 
employed in lattice calculations with the approach of per- 
turbation theory, which operates directly on the gauge- 
dependent elementary degrees of freedom. In addition, 
first hints how to provide such a relation have been ob- 
tained. 

This exploratory investigation has furthermore shown 
that the difference between the would-be confinement 
phase and the would-be Higgs phase is essentially man- 
ifest in the gauge-fixing sector outside the physical do- 
main, in agreement with previous studies and arguments 
p4l - |46| . It has also been shown that, within the limited 
range of gauges investigated here, there is little quali- 
tative difference between both phases otherwise. How- 
ever, the behavior found in the would-be Higgs phase 
is in agreement with the expectations for a dynamically 
induced Higgs effect. On the other hand, the confine- 
ment phase showed little difference compared to Yang- 
Mills theory. This is similar to the case when including 
dynamical fundamental quarks instead of scalars, where 
the difference haven been found to be also rather small 
(see, e.g., mm). 

In the present case, it is impossible to assess what 
the continuum and infinite- volume limit of these results 
will be, if the theory should be non-trivial. Therefore, it 
must now be the primary aim to further develop this ap- 
proach, and push in both phases towards the continuum 
field theory, or at least to a useful cutoff-theory in case 
of triviality. This will not be simple in this approach. 
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FIG. 5: The gauge dependence of the correlation functions in comparison between the Yang-Mills theory and the confinement 
phase. The top left panel shows the W propagator, the top right panel the ghost dressing function, the bottom-left panel 
the Higgs propagator (only the confinement phase is shown, see text), and the bottom right panel shows the running gauge 
coupling. The data points are in the minimal Landau gauge, and the hatched area represents the variability of the average 
value with the gauge. Diagonally hatched areas are for the confinement phase, and vertically hatched areas are for the case of 
Yang-Mills theory. For clarity, scale uncertainties have been suppressed in this figure, as they are not relevant since the ratio 
of scales of both cases is fixed. 



since the Higgs propagator itself will not be sufficient 
to differentiate between an only very weakly interacting 
theory and a trivial one. This requires determination of 
the vertices, a technique well-developed for Yang-Mills 
theory [24, 53], and straightforwardly applicable to the 
present setting. Also, results in the confinement phase 
are hopefully shedding light on how the confinement of 
fundamental charges in dynamical theories is manifest 
in low-order correlation functions [ir], and how it is 
different from screening due to the Higgs effect. 



Aclcnowledgments 

I am grateful to C. B. Lang and O. Rosten for help- 
ful discussions and to R. Alkofer, J. Greensite, and J. 
M. Pawlowski for comments on the manuscript. This 
work was supported by the FWF under grant number 
M1099-N16. The ROOT framework [sil has been used 
in this project. Computing time was provided by the 
HPC center at the Karl-Franzens-University Graz and 
the Slovak Grant Agency for Science, Grant VEGA No. 
2/6068/2006. 



10 



Effective t| mass 



'^ 4 

(0 



+ 

O 3 
O 



% 



- 










- 


• 








- 




T 
• 
1 




' 


" 
















- 








• Confinement phase 


" 


o 








- 




o 




o Higgs phase 








- 






§ < 


> 













t/a 



FIG. 6: The effective mass [5511 as a function of time. 



Confinement phase 

Fit 

IHiggs phiase 

Fit 
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the raw data together with a fit of the type (|A1[) . 



Appendix A: Determination of the scale and 
systematic uncertainty 

As noted in section lU the ground-state energy of the 
state (HI) is used to set the scale. This energy can be 
extracted from the exponential decay of its correlation 
function C{t) at asymptotically large times [531, '^^ P^^" 
ticular from a plateau in ln(C(i + 1)/C{t)). However, 
the correlator may also contain contributions from ex- 
cited states at finite times. As is visible in figure |6l this 
is indeed the case, at least in the Higgs phase. Thus, it is 
necessary to take the presence of excited states explicitly 



into account. A convenient way is to fit the correlator 
with the ansatz [53 

Cf{t) ^ acosh{m{t ~ N/2)) + bcoah{n{t - N/2)), (Al) 

where all quantities are in lattice units. The result, to- 
gether with the raw data, is shown in figure [T] Errors are 
determined by fitting besides the central value the one- 
a intervals, and as fit intervals always the time interval 
is used for which the correlation function is a positive, 
monotonously decaying function of time, and is adjusted 
for both upper and lower one-cr bounds. 

Despite the large number of configurations included, 
see table [TTl it is in particular for the confinement case 
not possible to unambiguously identify an excited state. 
Usage of an enlarged base of operators may here be use- 
ful in future calculations [55| . Therefore, only a single 
state fit is possible. The resulting fit parameters, to- 
gether with the number of configurations, are given in 
table [TTl The procedure in the main text is then used to 
associate a scale with the lower mass in the Higgs phase 
and the only mass in the confinement phase. The result- 
ing value is given in table |T| The error on the scale is just 
the statistical error. It is propagated through all calcula- 
tions in section [TTTI However, the dressing functions, and 
thus also the running coupling depending only on these, 
are dimensionless quantities, and therefore only the de- 
termination of the momentum in these cases receive an 
additional error. Furthermore, since the quenched case 
has a uniquely assigned scale based on the scale in the 
confinement case, the scale error is irrelevant when com- 
paring only the confinement and the quenched case, and 
can therefore be dropped in section [TVl 

Unfortunately, with the available lattice techniques it 
can never be guaranteed that there is not a further state 
still present, which is still lighter than the lightest one 
found, in particular if having a very small pre-factor. 
This is a source of systematic uncertainty in the scale 
determination. To illustrate this uncertainty, it will now 
be assumed that in the confinement phase there exists 
a lighter state with the same mass as the lighter state 
in the Higgs phase. For this hypothetical situation the 
error on this scale can be ignored. The results for the 
W sector in this case are shown in figure [51 and for the 
Higgs sector in figure [9l In the gauge sector nothing 
dramatic changes. Only the W propagator becomes less 
tree-level-like, while the running coupling becomes some- 
what more similar between both cases. Still, the ghost 
is the one most different between both cases, being a 
strongly momentum-dependent function in the confine- 
ment phase, but almost momentum-independent in the 
Higgs case. In the Higgs case, the changes are essen- 
tially negligible. Thus, the statements made in the main 
text are likely only (weakly) quantitatively affected by 
the systematic uncertainty in the scale determination. 
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FIG. 8: The W propagator (top panel), the ghost dressing 
function (middle panel), and the running coupling (bottom 
panel) under the assumption of equal scales for the Higgs and 
the confinement case. 
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TABLE II: The results for the fit of type (|A1[I . For the fits, the number of configurations indicated have been used, which have 
been generated in multiple runs with 1080 thermalization sweeps and 24 decorrelation sweeps. 
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FIG. 9: The Higgs propagator under the assumption of equal 
scales for the Higgs and the confinement case, the single large 
error bar is due to a numerical coincidence in the error propa- 
gation of the renormalization process, due to the appearance 
of a derivative. 
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